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ABSTRACT 

The relative contribution of baryons and dark matter to the inner regions of spiral 
galaxies provides critical clues to their formation and evolution, but it is generally 
difficult to determine. For spiral galaxies that are strong gravitational lenses, however, 
the combination of lensing and kinematic observations can be used to break the disk- 
halo degeneracy. In turn, such data constrain fundamental parameters such as i) the 
mass density profile slope and axis ratio of the dark matter halo, and by comparison 
with dark matter-only numerical simulations the modifications imposed by baryons; 
ii) the mass in stars and therefore the overall star formation efficiency, and the amount 
of feedback; iii) by comparison with stellar population synthesis models, the normal- 
ization of the stellar initial mass function. In this first paper of a series, we present 
a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for 
which multi-band high-resolution images and rotation curves were obtained using the 
Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on 
Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure 
systems. We characterize the sample of deflector galaxies in terms of their morpholo- 
gies, structural parameters, and stellar masses. We find that the SWELLS sample of 
secure lenses spans a broad range of morphologies (from lenticular to late- type spiral), 
spectral types (quantified by Ha emission), and bulge to total stellar mass ratio (0.22- 
0.85), while being limited to > lO^^-^ Mg. The SWELLS sample is thus well-suited 
for exploring the relationship between dark and luminous matter in a broad range of 
galaxies. We find that the deflector galaxies obey the same size-mass relation as that 
of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. 
We conclude that the SWELLS sample is consistent with being representative of the 
overall population of high-mass high-inclination disky galaxies. 

Key words: galaxies: spiral - galaxies: structure - galaxies: haloes - galaxies: fun- 
damental parameters - gravitational lensing 



1 INTRODUCTION 

The discovery of extended flat rotation curve s in the outer 
parts of disk gal axies tliree decades ago (|Bosmal 1 19781 : 
I Rubin et al.l 1 19781 ) was decisive in ushering the paradigm 
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shift that led to the now standard cosmological model dom- 
inated by cold dark matter (CDM) and dark energy. On 
cosmological (i.e. linear) scales there is excellent agreement 
between the predictions of the standard ACDM model and 
observations of the CMB, ty pe la Supernovae, w eak-lensing, 
and galaxy clustering (e.g. ISpergel et al.l [20071 ). On galac- 
tic and sub-galactic (i.e. non-linear) scales, however, there 
are a number of apparent discrepancies between the pre- 



2 Treu et al. 



dictions (e.g. iNavarro et al.lll997l: iBuUock et"aLll200lh and 



observations 
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2007 



Sand et al..,2008. : .Newton ct al.. .20111 ') of the structure and 



mass function of dark matter haloes. Whether these prob- 
lems reflect an incomplete understanding of galaxy forma- 
tion, or signal an unresolvable problem for the standard 
paradigm, remains to be determined. Measuring the den- 
sity profiles of dark matter halos on galaxy scales is thus a 
crucial test for the standard paradigm of galaxy formation, 
offering great potential for discovery. 

In most luminous galaxies the baryons (i.e. stars and 
cold gas) always make a non-negligible contribution to the 
gravitationa l potential, and due to the so-called disk-halo 
degeneracy (|van Albada fc Sancis il ll986l ) the distribution of 
baryons and dark matter in galaxies is poorly constrained, 
even when high spatial resolution rota t ion curves are avail- 
able (e.g. Ivan den Bosch fc Swatersl 1200 ll : iDutton et al.l 
l2005h . 

A potentially powerful method for breaking the disk- 
halo degeneracy is to co mbine galaxy kinem atics with strong 
gravitational lensing (|Maller et al.l |2000| ). This method 
makes use of the fact that kinematics measures mass en- 
closed within a sp here, while strong lensing measures pro- 
jected mass (e.g. iTreul [2010| . and references therein). The 
power of strong gravitational lensing and dynamics to break 
the disk-halo degeneracy has, to date, not been fully re- 
alized, due to the scarcity of known spiral galaxy grav- 
itational lenses. Up until a few years ago, only a hand- 
ful of spiral galaxy lenses with suitable inclinations to en- 
able rotation curve m easurements were known: Q2237+0305 
llHuchra et al.lll985l): B 1600+434 (|jaunsen fc Hiorthlll997l: 
Koopmans et al.1 1 19981 ): PMN J2 004-1349 (IWinn et all 



2003h : CXOCY J220132. 8-320144 (|Castander et al.1 120061 ). 



However most of these systems are doubly imaged QSOs, 
which provide minimal constraints on the projected mass 
distribution. Q2237H-0305 is a quadruply imaged QSO, 
which gives more robust constraints, but since the Ein- 
stein radius is five times smaller than the half light size 
of the galaxy, the lens ing is mostly sensitive to the bulge 
mass iTrott et al.ll201ol ). Recently, extensive searches based 
on large parent surveys have started to find gravitational 
lens candidates where the deflector is a spiral galaxy. Ex- 
amples of imag ing-based searches f or spiral lens galaxies in - 
clud e those by ISygnet et al.1 (|201G| ) . [Marshall et al.1 (|2009t ). 
and iMore et al.l mill), while spect roscopy-based searches 



incl ude those by 
vey (|Bolton et al. 



Feron et al. 



20061. 2001 



(I2OO9I) and the SLAGS Sur- 



Auger et al]|2009l ). One of the 



main challenges of identifying and studying strong lens sys- 
tems where the deflector is a spiral galaxy is the presence 
of dust. Modeling and correcting for dust obscuration is es- 
sential to model the gravitational potential of the deflector; 
this can be done by exploiting the conservation of surface 
brightness of the background source in gravitational lensing. 

Since its beginning in Ifubble Space Telescope (HST) 
Cycle 13, SLAGS has discovered 85 strong galaxy-galaxy 
lenses (and an additional 13 high probability candidates). 
Because of the strong dependency of the lensing cross- 
section on central surface mass density, ~ 90% of the lenses 
disc overed by SLAGS have been massive early-type galax- 
ies (|Auger et al.ll2009l ). However, thanks to its large overall 
size, SLAGS was able to discover a considerable number of 



previously unknown gravitational lens systems with a de- 
flector of spiral morphology. 10 SLAGS lenses have visible 
spiral morphology (classification S or SO/Sa). Of these 10, 8 
have disk inclinations high enough {b/a < 0.6) for measur- 
ing reliable rotation curves. This subset alone already more 
than doubles the sample of previously known systems. In 
addition, source and deflector redshifts are known by con- 
struction for the SLAGS systems, and the deflector galaxy is 
generally bright and at a relatively low redshift, making it a 
very practical sample for detailed follow-up. The main limi- 
tations of the SLAGS-spiral sample are that it is limited to 
bulge-dominated spirals, since they are the most common in 
strong lensing selected samples, and that spatially resolved 
kinematics is not available from the SPSS fiber-based spec- 
troscopy, but requires follow-up (|Barnabe et al.|[201ll ). 

The goal of the SWELLS (Sloan WFG Edge-on Late- 
type Lens Survey; WFG = Wide Field Gamera) project is 
to identify, follow-up, and analyze in a systematic fashion 
a large and homogeneous sample of spiral lens galaxies. We 
aim to combine strong gravitational lensing with stellar kine- 
matics to break the "disk-halo" degeneracy and address two 
fundamental issues. Firstly, measurements of dark matter 
halo density profiles and halo shapes will be compared to 
predictions of galaxy formation models in the concordance 
AGDM cosmology. Secondly, measurements of stellar mass- 
to-light ratios will be compared to those predicted from stel- 
lar population models, t hus placing constraints on the stel- 
lar initial mass funct ion (|Treu et al.l2010l : lAuger et al bold : 
ISpiniello et"al]|201ll ). 

The sample consists of high-inclination spiral lens galax- 
ies, drawn from both the SLAGS Survey and a new dedi- 
cated SDSS/i/5T search. For each object in its target list, 
the SWELLS survey is measuring high quality rotation and 
velocity dispersion profiles using long-slit spectroscopy at 
the Keck Telescopes. In addition, high resolution images in 
multiple bands, including the infrared with Keck, are being 
obtained to model and correct dust obscuration, identify 
and model multiple images of the background source, and 
infer stellar masses of bulge and disk from stellar population 
synthesis models. 

In this first paper of the series we present an overview 
of the survey and the lens sample, and discuss the lensing 
galaxies as spiral galaxies to identify potential selection ef- 
fects that must be taken into account when interpreting our 
results. Papers of this series will present detailed analyses 
of the individual systems illustrating the methodology and 
challenge s of a joint lensing and dynamical study of spiral 
galaxies (|Dutton et al.l 1201 ll . hereafter paper II) and a li- 
brary of gravitational lens models (Brewer et al. 2011, in 
preparation) . 

This paper is organized as follows. Our strategy and 
sample selection algorithm is described in Section (2] The 
data are described in Section (3] In Section |4] we present 
lensing classification, morphological and structural proper- 
ties and stellar masses for each system. In Section [5] we 
discuss the lensing galaxies as spiral galaxies, by compar- 
ing them to a control sample drawn from the SDSS survey. 
We discuss our findings and the lessons learned with an eye 
to future searches for strong lensing spirals in Section |6] A 
brief summary is given in Section [T] Throughout this paper, 
and the rest of t he SWEL LS series, magnitudes are given in 
the AB system (|Okelll973 ) and we adopt standard "concor- 
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Table 1. Summary of basic properties of the SWELLS sample of spiral lens galaxies. Col. 1 lists the lens ID; Cols. 2 and 3 the 
coordinates (J2000); Cols 4 and 5 give deflector and sour ce redshifts; Col. 6 the lens grade (A=socure, B=probable, C=possible, X=not 
a lens) according to the definition of l lBolton et al.|[2"00(i) : Col. 7 the velocity dispersion within the spectroscopic fiber as measured by 
SDSS; Col. 8 the Ho flux as measured by SDSS spectroscopy based on the subtracted model spectrum; Col 9. the axis ratio as measured 
by SDSS; Col. 10 the reference for the discovery of the lens (l=Bolton et al. 2006, 2=Bolton et al. 2008, 3=this paper). 



dance" cosmological parameters, i.e. h = 0.7, Q,m. = 0.3 and 
Q.A ~ 0.7, where the symbols indicate the Hubble Constant 
in units of 100 km s~^ Mpc~^ and the matter and dark en- 
ergy density of the Univer se in units of the critical density 
(e.g. iKomatsu et"al]|2009l ). 



2 SAMPLE SELECTION 

The SWELLS sample is composed of two subsamples. The 
first is a subset of the spiral lens galaxi es identified by th e 
SLAGS Survey (|Bohon et al. 2006, 2008; Auger et al]|2009l ). 
suitable for kinematic follow-up due to their high inclination, 
as determined by the observed axis ratio. The second sub- 
sample consists of new spiral lens galaxies, identified from 
the SDSS-database and confirmed with HST imaging in the 
following manner. The strategy is similar to that adopted 
for the original SLAGS sample, optimised to identify spiral 
lens galaxies. 

First we looked for SDSS spectra that contain two sets 
of lines at two different redshifts, indicating a foreground 
and background object in the same 3 arcsecond-diameter 
solid angle covered by the fiber as detailed in the Appendix. 
Out of almost a million spectra, our spectroscopic search 
algorithm found more than 200 new high-probability lens 
candidates. Second, SDSS images of the lens candidates were 
visually inspected by two of us (TT and AAD) to identify 
those spiral galaxies sufficiently inclined for rotation curve 
measurements, and to reject obvious failures or contami- 
nants. Third, we estimated the strong lensing probability 
according to the following procedure. Based on our previous 



experience ([Bolton et al ] |2008l ). the lens confirmation rate is 
a well known and monotonically increasing function of Ein- 
stein radius (effectively, proportional to the fraction of the 
solid angle of the fiber that is contained within the critical 
lines, plus seeing and other observational effects). The Ein- 
stein radius was estimated from the velocity dispersion of 
the lens and the usual ratio of angular diameter distances, 
assuming that the stellar velocity dispersion equals the nor- 
malisation of the best fitting singu lar isothermal ellipsoid 
(|Treu et al.ll2006l : iBolton et aLll200l ). Since velocity disper- 
sions are not always available or reliable, especially for the 
more disk-dominated systems with emission lines, we used 
total luminosity as a proxy for stellar velocity dispersion, af- 
ter calibrating the correlation on the rest of the sample. We 
expect this to be a conservative estimate, since edge-on sys- 
tems are mor e effective lenses due to their higher projected 
mass density (iMaller et al.||l997l : lKeeton fc Kochanekill998l : 
iBlain. Moller. fc Mailer 19991 '). Fourth, to optimise use of 
telescope time, we excluded all targets with estimated lens- 
ing probability below 30%. This left a total sample of 43 can- 
didates, which we used to select targets for high-resolution 
imaging follow-up based on observability. If more than one 
candidate was observable at a given time, we gave priority to 
those with the highest estimated lensing probability. Follow- 
up observations and lens confirmation rates are discussed in 
the next sections. 



A total of 27 candidates were observed with high resolu- 
tion imaging. A summary of the properties of these SWELLS 
lenses and candidates is given in Table [T] 
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3 OBSERVATIONS 

3.1 Hubble Space Telescope imaging 

In the supplementary Cycle 16s we were awarded 91 or- 
bits (of which 74 were executed) of HST time to discover 
new spiral galaxy lenses and to complete three-band opti- 
cal imaging of the known SLACS spiral lenses (GO 11978; 
PI Treu). As a compromise between maximum wavelength 
coverage - needed to model and correct dust obscuration, 
identify strongly lensed features, and model stellar popula- 
tions - and efficiency, we chose to image each lens through 
the F450W, F606W, and F814W fihers. 

As part of GO-11978, twenty new candidates (se- 
lected as described in the previous section) were sched- 
uled to be observed. Of those, fifteen were observed with 
WFPC2 in 3 bands: F450W (two orbits), F606W (one 
orbit), and F814W (one orbit), and one was observed 
in F450W and F606W (two and one orbit, respectively). 
Due to scheduling constraints the F814W exposure and 
remaining 4 targets were not observetQ- In addition, five 
SLACS-selected spiral lenses (J0841+3824, J1029-f0420, 
J1032-H5322, J1103-H5322, J2141-0001) were imaged in 
F450W (two orbits each) and/or F606W (one orbit) to com- 
plete three-band imaging for the sample. 

All targets were placed in the centre of the WFPC2- 
WF3 chip, which has a pixel scale of O'.'l. The observations 
in each filter consisted of four exposures dithered by frac- 
tional pixel amounts to remove defects, identify cosmic ray 
hits, and to recover resolution lost to undersampling. The 
raw images were reduced and combined using a custom built 
pipeline based on the drizzle package, resulting in final im- 
ages with a pixel scale of O'.'OS. 

Additional HST imaging is available for the SLACS sub- 
sample as listed in Table [U ACS/WFCP2 a nd NICMOS 
observations are presented in the p apers by iBolton et al.l 
l|2006l. l2008h and I Auger et all (|2009l ). Infrared images ob- 
tained with WFC3 through filter F160W as part of GO- 
11202 (PI Koopmans) are presented here for the first time. 
As with the other instruments, a standard four-point dither 
strategy was adopted, in this case to remove defects and 
recover some of the information lost to the undersampling 
of the point spread function (PSF) by the 0'.'13 pixels. The 
data were reduced via a custom pipeline based on the drizzle 
package, using the most up to date calibration files. 

3.2 Keck Adaptive Optics Imaging 

High-resolution infrared imaging is extremely valuable for 
dust correction, multiple image identification and stellar 
mass estimation. For these purposes, K'-hanA (2.15/im) 
imaging was obtained for a subset of candidates/lenses for 
which a suitable bright nearby star is available for tip- 
tilt correction using the Laser Guide Star Adaptive Optics 
(LGSAO) system on the Keck-II 10m Telescope. Higher pri- 
ority was given to those systems for which HST infrared 
imaging is not available, although the combination of HST- 
F160W and ground based if'-band imaging is especially 

^ Three of these remaining targets were observed with K-band 
laser guide star adaptive optics (LGS-AO) imaging on the Keck 
II telescope as described in § 13.21 



powerful. The sensitivity of Keck AO at longer wavelengths 
complements very nicely the HST PSF stability and sensitiv- 
ity to low surface brightness features, helping in particular 
with dust correction and stellar mass estimates. 

All observations were performed with the NIRC2 cam- 
era in wide field mode, yielding a pixel scale of O'.'Ol, typical 
resolution of O'.'lO FWHM and strehl ratio of 0.15-0.20. The 
observations took place on the following dates, during times 
of clear sky and good seeing conditions suitable for AO cor- 
rection: Nov 11 2009 (PI Koo); May 7/8 2009, August 12 
2009, January 04 2010 (PI Treu). Typical exposure times 
ranged between 30 and 60 minutes, depending on condi- 
tions and faintness of the target. Table [5] lists the targets 
that were observed with AO. More details on the observa- 
tions of individual systems will be given in future papers 
along with the analysis of the systems themselves. Although 
in general the lensed features are rather blue and more easily 
identified in HST images, we note that two of the grade-A 
SWELLS lenses were observed only with Keck- AO, and con- 
firmed from the ground. 

All NIRC2 data were reduced with a custom Python- 
based reduction package described bv lAuger et al.l (|201ll ). A 
sky frame and a sky flat were created from the individual sci- 
ence exposures after masking out all objects. Frames were 
then flat-fielded and sky-subtracted. The images were de- 
warped to correct for known camera distortion. The frames 
were aligned by centroiding on objects in the field; each 
frame was then drizzled to a common output frame, and 
these output frames were median-combined to produce the 
final image. 

An empirical model for the PSF was derived from ob- 
servations of stars. In some cases a star was close enough to 
the lens galaxy to fall in the NIRC2 field of view and be used 
as a model PSF. In other cases a PSF star pair was used; 
this separate observation involves observing a star cluster, 
and using as a PSF model a star that is at the same dis- 
tance from its tip-tilt correction star as the lens galaxy was 
from its tip tilt star. The star pair observations were made 
immediately following the lens observations. 



4 SAMPLE PROPERTIES 

In this section, we describe the overall properties of the sam- 
ple. First, in Section r4.1l we present our classification of the 
lensing morphology, providing notes on each of the systems 
based on visual inspection of multi-colour images and deflec- 
tor subtracted images. These images are presented in Fig- 
ures[l]and[2l We note that subtracting the light of the deflec- 
tor galaxy is much more challenging for the spiral galaxies 
of the SWELLS sample, compared, e.g., to the smoother 
early-type galaxies that compose the SLACS sample. For 
this reason the deflector subtracted images shown in this 
paper are used only for the classification of lensing morphol- 
ogy. In future papers, we will discuss alternative strategies 
to disentangle light from the source and the deflector for the 
purpose of lens modeling. 

Second, in Section 14.21 we describe the morphological 
and structural properties of each deflector galaxy, based on 
visual inspection and of multi-component, multi-band sur- 
face photometry. Third, in Section l4^ we derive stellar mass 
for the bulge and disk component of each of the deflectors. 
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We conclude in Section |4]4] by investigating the dependency 
of the lens confirmation rate on the properties of the deflec- 
tor and the lensing geometry. 

4.1 Lensing classification 

A summary of the SWELLS sample is given in Table [T] All 
lenses were given a classification based on the scherne and 
criter ia defined by the SLAGS survey (|Bolton et al.l |2006| . 
l2008l : I Auger et all 120091 ). By consensus of a subset of the 
authors (TT, AAD, MWA, PJM, BJB), lenses were classi- 
fied as secure (A-grade), probable (B-grade), possible (C- 
grade) and not a lens (X-grade). We note that our classi- 
fication scheme is quite strict, and requires the convincing 
identification of multiple images. It is likely that most of the 
B-grade lenses and some of the C-grade systems could turn 
into A-grade systems if deeper, or higher resolution, data 
were available to identify multiple images. 

In addition to the known late-type high-elongation 
lenses from the SLAGS survey (pre-Cycle 16s), 8 A-grade 
previously unknown strong gravitational lenses were discov- 
ered as part of this project. Of the remaining eleven sys- 
tems, there is 1 B-grade and 6 C-grade systems. 4 systems 
were classified as definitely not-lenses. Full gravitational lens 
models for the secure systems will be presented in future 
papers in this series (paper II; Brewer et al. 2011, in prepa- 
ration). 

4.1.1 Notes on individual systems 

0820+4847 Grade A. Three component arc to NE of bulge. 
Defiector has almost edge-on disk. Bulge is small, suggesting 
a later type spiral, although spiral arms are not visible due 
to inclination. Observed with LGSAO, source appears as two 
lensed images, possibly indicating that the source has a red 
bulge and an extended blue disk. 

0821+1025 Grade G. Observed only with LGSAO. 
Foreground galaxy is high inclination disk. Background 
galaxy is compact, 1.7" to the W. The background galaxy 
might be strongly lensed at fainter surface brightness levels. 
Blue band imaging might reveal multiple images. There is 
also another source, possibly a satellite of the main spiral 
galaxy, ~ 3" SW of center. 

0841+3824 Grade A. The deflector is a grand design 
spiral with blue arms and a redder bar/disk. A simple lens 
model, reproducing the curved lensed feature, ~ 1.5" NW of 
center and the counter image on th e opposite side, is given 
in the paper bv lBolton et all (|2008l '). 

0915+4211 Grade A. Observed only with LGSAO. 
Glearly visible arc and counter arc. 

0955+0101 Grade A. The deflector is an edge-on disky 
galaxy. A blue and highly sheared arc is visible to the SW 
of the bulge and the counter-image is clearly visible in the 
residuals, especi ally the B-band. A s imple lens model is given 
in the paper bv iBolton" et all (|2008l 'l. Another source, possi- 
bly a satellite of the defiector, or of the background source, 
is visible ~ 5" to the S of the defiector. 

1015+1750 Grade X. The defiector is a dusty star 
forming edge-on spiral. The likely background galaxy is 1.6" 
to SE. No F814W image. No evidence for highly sheared im- 
ages. 



1029+0420 Grade A. The deflector is smooth and 
disky and seen edge-on, likely an early-type spiral or lentic- 
ular galaxy for the relatively large bulge. Arc and counter 
arc are clearly visible in the colour image and in the bspline 
su btracted imag e . A si mple lens model is given in the paper 
bv lBolton et all (|2008l ). 

1032+5322 Grade A. The deflector is smooth and 
disky and seen edge-on. The lensed source is blue and has 
a complex morphology, likely a spiral with a redder bulge. 
Images and counter- images are clearly visible. A sim ple lens 
model is given in the paper bv lBolton et al.l ()2008l ). 

1037+3517 Grade A. The main deflector is a spiral 
galaxy. However, a second peak in red surface brightness in- 
dicates either a merger in process or a multi-galaxy lens, if 
the two galaxies are not in the same plane. Arc and counter 
arc are blue and clearly visible, although the morphology of 
the lensed images is not straightforward to interpret, owing 
to the complex lens potential. Gould be a multi-galaxy lens, 
possible 2nd (compact red) lens galaxy at 1.1" to N. Ex- 
tended source galaxy 2.4" to N. Integral field spectroscopy 
or high resolution narrow band imaging would help decipher 
this configuration. 

1056+0005 Grade X. Defiector galaxy is a bulge-less 
star forming spiral. Likely background galaxy is seen l"3 
North of centre. No evidence of shear suggests that the 
source is well outside the caustics of the deflector. 

1103+5322 Grade A. The deflector is smooth and 
disky and seen edge-on. The lensed source is blue and has 
a complex morphology, likely a spiral with a redder bulge. 
Images and counter-images are clearly visible. A sim ple lens 
model is given in the paper bv lBolton et all (120081 ). 

1117+4704 Grade A. The deflector is smooth and 
disky and seen close to edge-on. A bright and highly sheared 
source is visible l'.'6 NE of deflector centre, albeit with no 
clear evidence for a counter-image. A classic "fold" quadru- 
ply imaged source is visible in the B-band and V-band b- 
spline residuals, with radius approximately 0'.'75. 

1135+3720 Grade A. Observed only with LGSAO. Arc 
and counter-image are clearly visible in the b-spline residu- 
als. 

1203+2535 Grade A. Deflector galaxy has a blue edge- 
on disk and a red bulge; the outer edge of the disk is warped. 
Arc made of three merging images is clearly visible SE of 
bulge. Possible quad with faint counter-image near the cen- 
tre or naked cusp configuration. 

1228+3743 Grade B. Defiector is a grand design spiral. 
Background source is also extended with a red bulge and 
blue spiral arms, one of which passes close to the lens bulge. 
A possible arc and counter arc with image separation ~ 2" 
is visible in the residuals, clearest in the B-band. 

1251-0208 Grade A. The defiector has a red bulge and 
blue spirals arms. The lensed source is blue and high sur- 
face brightness. Images and counter-images are clearly visi- 
ble in the b-sp line residuals . Simp l e lens models are given in 
the papers bv lBohon et al.l (|2006l . l2008l ). A l ensing and dy- 
namic s analysis is presented in the paper bv lBarnabe et ah! 
(|201lh . 

1258-0259 Grade X. Foreground galaxy has edge-on 
disk and bar like structure. A possible background galaxy is 
seen l'.'6 S of the defiector, with bulge and spiral arms. The 
source does not appear to be sheared and no counter-image 
candidate could be identified. 
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Figure 1. Colour montage of the high resolution images of all the SWELLS targets. The image sizes vary from 10" to 1' on a side, 
and each lens' grade is shown in parentheses next to its name. The three systems with grey-scale images only have AO K-band imaging. 
North is up, East is left. 
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Figure 2. Single band high resolution images of the SWELLS targets. For each system the top row shows the original image, and 
the bottom row shows the residual image after B-spline subtraction of the deflector. Missing bands are replaced by a cross. The top 
images are 10" on a side, while the residual images are 5" on a side to show the details of the potential sources identified in the SDSS 
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1300+3704 Grade X. Foreground galaxy is a large 
bulge-less edge-on star forming spiral. Possible background 
galaxies 3'.'5 to NE and SW of deflector. 

1313+0506 Grade A. The deflector is smooth and 
disky and seen edge-on. The lensed source is blue and high 
surface brightness. Images and counter-images are clearly 
visible in the residuals . A sim ple lens model is given in the 
paper bv lBolton et all (|2008l '). 

1321-0115 Grade C. Foreground galaxy has edge-on 
disk with prominent dust lane. Blue background galaxy is 
seen 1'.'9 to N of the deflector. The source is elongated tan- 
gentially to the direction of the deflector suggesting strong 
shear. No counter-image could be identified, although it 
could be obscured by the prominent dust lane. High res- 
olution infrared images would be useful to confirm or reject 
the lensing hypothesis. 

1331+3638 Grade A. Deflector galaxy has blue spiral 
arms and a red bulge. Arc made of three merging images is 
clearly visible NE of bulge. Possible quad with faint counter- 
image near the centre or naked cusp configuration. The tar- 
get was observed with LGSAO, but faint tip-tilt star yielded 
poor correction. 

1410+0205 Grade C. The deflector has a red bulge, 
a prominent disk and a ring-like blue structure surround- 
ing the bulge. The blue structure is most likely due to star 
formation in the deflector although lensing cannot be ruled 
out with current data. Another possible source is seen NW 
of lens centre, although direction of the elongation is not ob- 
viously consistent with strong lensing. It could be that the 
source is only partially multiply imaged with the bulk of it 
extending outside of the caustic. Integral field spectroscopy 
or narrow band imaging would be useful to decipher this 
system. 

1521+5805 Grade C. Foreground galaxy has nearly 
edge-on disk with prominent dust lane. The background 
galaxy, with bulge and spiral arms, is visible NW of de- 
flector. It is possible that parts of the background source 
are multiply imaged but the dust lane and the complex 
morphology prevent a conclusive answer. High resolution in- 
frared images and integral field spectroscopy or narrow band 
imaging would help. 

1556+3446 Grade C. The deflector is a dusty nearly 
edge-on spiral galaxy. A possible image and counter-image is 
visible in the b-spline residuals. In the SDSS spectrum, back- 
ground source has clear Oil, Hb, OIII emission lines. High 
resolution K-band imaging would help see through the dust, 
while narrow band imaging or integral field spectroscopy 
would help disentangle the light from the source and from 
the deflector. 

1703+2451 Grade A. The deflector is an edge-on dusty 
spiral, with boxy isophotes suggesting a pseudo bulge. The 
outer disk is warped. A red(dened) arc and possible counter- 
image are clearly visible, especially in the LGSAO K-band 
b-spline residuals. 

2141-0001 Grade A. Described in det ail in paper 
II (iDutton et al] 1201 ll ) , and more simply by iBolton et al.l 
(l2008l^ . 

2210-0934 Grade C. The deflector has a prominent 
edge-on disk. Possible image and counter-images are visible 
in the b-spline residuals, with image separation ^ I'.'S. An- 
other source, possibly a satellite of the main lens, is visible 
~ 5" N of center. 



4.2 Morphological and structural properties of 
the deflectors 

All the SWELLS deflectors show a clear disky component 
(Figure [T|, confirming our original selection based on SDSS 
images. Furthermore, the targets have high projected ellip- 
ticity, again consistent with our original selection based on 
SDSS images (Figure|3}. Many of the targets show clear spi- 
ral arms, which in several cases are blue. Other targets are 
so close to edge-on that it is hard to identify spiral arms 
conclusively (notes on the morphology of individual targets 
are provided in the previous section). The overall morphol- 
ogy of the targets is also consistent with the distribution of 
velocity dispersions (Figure this distribution has a lower 
central value than the SLAGS sample, which is dominated 
by massive elliptical galaxies. 

We quantify the structural properties of the deflectors 
by performing a bulge-disk decomposition. For the purpose 
of this paper we adopt a simpl e parametrisation, an d de- 
scribe the bulge with an elliptical lde VaucouleursI (|l948l ) pro- 
file and the disk with an elliptical exponential profile. This 
parametrisation is the same as that adopted by Simard et 
al. (2011, in preparation) to analyse the entire SDSS-DR7 
database, and therefore is ideal for a comparison between 
lenses and non-lens galaxies. Beyond this first analysis, in 
future papers we will derive more accurate descriptions of 
the surface brightness distribution by considering more flex- 
ible models, like, e.g., the Sersic profile for the bulge (e.g. 
paper II) , and by correcting for extinction due to prominent 
dust lanes (Brewer et al. 2011, in preparation). 

For each deflector we compute the posterior probability 
distribution function of the parameters by comparing the 
models with the multi-colour images using the dedicated 
code SPASMOID d eveloped by one of us (MWA; see also 
iBennert et"al]|2011^ . SPASMOID replaces the functional- 
ity of standard codes like GALFIT, by enabling a full explo- 
ration of the likelihood surface as well as joint analysis of 
multiband data. Furthermore, SPASMOID models fit the 
multi-band imaging simultaneously with the same effective 
radius, and the bulge effective radius is constrained to be 
smaller than the disk effective radiufl This procedure en- 
sures self-consistent colors even with data of varying depth 
as a function of wavelength. The most probable values of 
the parameters are given in Table [2] 



4.3 Stellar Mass 

Stellar mass estimates for the bulge and disk compo- 
nents of each candidate were derived by comparing multi- 
band photometry w ith stellar population synthesis models 
j Auger et al.ll2009l ). The method we employ provides not 



^ The priors for the surface brightness models are uniform but 
have meaningful bounds: ±5 pixels from the brightest central 
pixel for the centroid, 0.1 (0.15) to 1 for the axis ratio of the 
disk (bulge), and effective radii between 0.1 arcsec (2 pixels) and 
the size of the cutout that is fitted by SPASMOID (typically hun- 
dreds of pixels). Furthermore, the effective radii, flattening, and 
position angle are the same for each filter, while the centroids 
of the bulge and disk are allowed a small shift of 3 pixels with 
respect to each other. 
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Table 2. Photometric, structural and stellar mass parameters of the SWELLS systems. For each lens the first line lists the bulge 
parameters (when detected), while the second line lists the disk parameters. Ty pical uncertainties are 0.05 mags on magnitudes, 10% 
on effective radii and a^xis ratios (q = b/a), dominated by systematics (see, e.g.. I Auger et al]|2009l , and references therein). Apparent 
magnitudes are not corrected for galactic extinction. Effective radii are given along the intermediate axis; major-axis effective radii can 
be obtained by dividing by 



only stellar masses but also uncertainties that take into ac- 
count all the relevant degeneracies between, e.g., age and 
metallicity of the stellar populations. The dominant system- 
atic uncertainty is the unknown normalisation of the stellar 
initial mass function (IMF). Determining the absolute nor- 
malisation of the IMF is one of the goals of the SWELLS 
survey (see paper II). For the moment, we adopt two pop- 
ular choices for the IMF, which are believed to bracket th e 
appropriate range for s piral galaxies |Bell fc de Jons)[200ll) 



-^^--^ „ s piL__ - - - ■- ,.- 

Stellar masses for these lchabrle^ (|2003l ') and lSalpeteJ (|l955l ) 
IMFs are given in Table [3 



4.4 Notes on the selection function 

The overall success rate for confirmation of SWELLS can- 
didates in Cycle 16s is 42 ± 13% (8/19) for A-grade and 
47 ± 14% (9/19) for A-|-B-grade lenses. As expected, this 
is somewhat lower than the confirmation rate for massive 
early-type galaxies in the SLAGS survey, due to the fact 
that our galaxies are on average less massive: the Einstein 
radii are smaller and occupy a smaller fraction of the SDSS 
fiber. Furthermore, identifying multiple images is inherently 
more difficult behind spiral deflectors, due to the complexity 
of their surface brightness distribution and the presence of 
dust. 
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Figure 3. Histogram of SDSS parameters for SWELLS lenses. 
Tile top panel sliows tlio distribution of stellar velocity disper- 
sions (when available, see Table [TJ, the middle panel shows the 
distribution of deflector redshifts, and the bottom panel shows 
the distribution of axis ratios. Systems discovered prior to Cy- 
cle 16s are plotted in red, our new systems are shown in black. 
Solid histograms represent secure lenses (A), hatched histograms 
represent probable and possible lenses (B and C), while open his- 
tograms represent non-lenses (X). Note that our Cycle 16s strat- 
egy targeted preferentially lower velocity dispersion and axis ratio 
systems, in order to select disk edge-on galaxies. Note also how 
the confirmation rate increases with cr, as expected. 



We can investigate the selection function more quanti- 
tatively by studying the stellar mass structure of the deflec- 
tors and correlations of its properties with the confirmation 
rate. Mass and surface mass density are expected to be the 
key parameters determining the lensing strength of a galaxy. 
We will therefore investigate the dependency of confirmation 
rate on several proxies for mass and density. 

First, in Figure |4] we show the average stellar mass 
within circular apertures, expressed in units of the critical 
density Ecrit = 4^gji j — i where Ds and Dds are the an- 
gular diameter distance to the source, to the deflector and 
between the deflector and the sou rces, respectively (see, e.g., 
ISchneider et al.ll 19921 : iTreulbOld ). As shown in Figure [1 the 
enclosed average stellar mass density is always dominated 
by the bulge in the inner parts, although the disk starts to 
contribute significantly beyond 1 kpc. All of the secure sys- 
tems have deflectors of uniformly high density, above critical 
density well beyond one kpc, even for a Chabrier IMF. The 
not-lenses span a broader range in surface mass densities 
including several systems that are significantly underdense 
compared to the grade-A systems. Three out of four of the 
X-grade lenses do not have a detectable bulge component 
(to our sensitivity limit of M, = lO^M©). 

In a circularly symmetric system the Einstein radius 
corresponds to the radius within which the average m ass 
density equals the critical density (|Schneider et aLlliggd ). In 
order to obtain a proxy for the Einstein radius even when 
no lensing is detected we define the "stellar Einstein radius" 
as the radius within which the average stellar mass density 
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Figure 4. Average (enclosed) stellar surface mass density profile 
of SWELLS targets, in units of the critical density. The "stellar 
Einstein radius" of each lens lies at the point where (S)/E(,i-it — 
1. The black lines represent the total mass density profile, the 
red lines represent the contribution of the bulge, and the blue 
lines represent that of the disk. Profiles for both the Salpeter 
and Chabrier IMF are shown. The top panel shows the profiles of 
secure (A) lenses, the middle panel shows the profile of probable 
(B) and possible (C) lenses, and the bottom panel shows the non- 
lenses (X). Note that the bulge always dominates the enclosed 
average mass at small radii, while the disks start to orovide a 
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B/T (M.) 

Figure 5. Distribution of stellar mass related properties. The top 
panel shows the distribution of stellar Einstein radii, i.e. the radii 
within which the stellar surface mass density is equal to 1 in criti- 
cal units, assuming a Chabrier IMF. This is an approximate lower 
limit to the actual Einstein radius, which will include the contri- 
bution of dark matter. The middle panel shows the distribution 
of total stellar mass. The bottom panel shows the distribution 
of bulge to total stellar mass ratio (the "bulge fraction"). As in 
Figure [S] the red histogram represents lenses discovered prior to 
Cycle 16s; the black solid, hatched and empty histograms repre- 
sent lenses targeted in Cycle 16, identified respectively as secure, 
probable and possible, and not-lenses. Note how no lens was con- 
firmed with stellar Einstein radius below 0"94. 



equals the critical density. Note that high resolution multi- 
band imaging is necessary to estimate this quantity robustly. 
Clearly this is only an approximation for non-circular sys- 
tems, and will be a lower limit to the true Einstein radius if 
dark matter is present. 

We are now in a position to study the distribution of 
integrated stellar mass M*, bulge-to-total stellar mass ra- 
tio B/T and stellar Einstein radii for the SWELLS targets. 
As shown in Figure [5] all three quantities correlate at some 
level with confirmation rate. Most of the not-lenses are found 
at small M, , small B/T, and small stellar Einstein radius. 
However, the correlation is not perfect: B and C grade lenses 
are found even at M* > 10^^ and B/T> 0.7. All three ob- 
jects with no detected bulge turned out to be not-lenses. 
Perhaps the cleanest predictor of confirmation is stellar Ein- 
stein radius: there are no confirmed lenses below 0'.'94 even 
though there are probable and possible lenses well above 1". 
For comparison, the stellar Einstein radii for the confirmed 
SLAGS early-type deflector galaxies are also typically (90%) 
above 0'.'94 although 10% of the values are in the range 0'.'69- 
0'.'94. This tail of the distribution extending to smaller radii 
is consistent with the idea that multiple images with smaller 
separation can be detected in the presence of smooth early- 
type deflector galaxies compared to the case of spiral deflec- 
tors. However, the number of SWELLS targets is too small 
too reach a statistically significant conclusion. 

In conclusion, the selection function of the new 
SWELLS targets is complex and comprised of multiple 



steps. In the first one, we preselect based on SDSS images 
edge-on late type galaxies. This is extremely efficient in pro- 
ducing a sample of lenticular and spiral galaxies. The second 
step, confirmation as a lens system, depends on the prop- 
erties of the deflector but also has a stochastic component. 
Denser systems produce larger Einstein radii, which in turn 
make the lensing effect more likely and easier to detect at 
HST resolution. In addition, the relative position on the sky 
of the source and deflector, as well as the broad band surface 
brightness of the source, ultimately control the presence of 
detectable multiple imaging. A third layer of complication 
is added by the presence of dust lanes and small scale struc- 
ture in the surface brightness distribution of the deflector, 
which may complicate identification of multiple images even 
for a perfectly good strong lensing system. In the next sec- 
tion we will investigate the effective selection function of 
the SWELLS sample by comparing the properties of the 
SWELLS deflectors with those of a typical sample of spiral 
galaxies. 



5 SWELLS GALAXIES AS SPIRAL GALAXIES 

Understanding the selection function of any observational 
sample is essential, both in order to generalise results to its 
parent population, as well as to compare across samples. In 
the case of the SWELLS sample the selection function of the 
deflectors is sufficiently complicated that it is extremely hard 
to compute it accurately from first principles. We therefore 
take an alternative approach and reconstruct the basic fea- 
tures of the selection function a posteriori by comparing the 
properties of the SLAGS detectors to those of a suitable com- 
parison sample of non-lens galaxies. The comparison sample 
should be as close as possible in terms of controllable galaxy 
parameters, so that significant differences might be inter- 
preted unambiguously as due to the lensing selection effect. 
It is particularly useful to perform this comparison for both 
confirmed lenses and not-confirmed lenses as this will help 
separate the effects of our SDSS-based pre-selection from 
those related to the ability of the deflectors to produce de- 
tectable strong lensing events. In this paper we focus on the 
properties of the SWELLS deflector galaxies, and leave for 
future work the issue of whether they live in an overdense 
environment (whi ch turns out not to be the c ase for the 
SLAGS lenses, see I Auget|[200^ : iTreu et al.ll2009l V 

Among the quantities that are well-studied and avail- 
able for large samples of galaxies, stellar mass and size (and 
therefore density) are particularly interesting and likely to 
influence the lensing selection function. We therefore choose 
to focus on the size-stellar mass correlation as our compari- 
son tool (for a similar analysis of the SLAGS sample, b ased 
on the fundamental plane relation see iTreu et al.ll2006l '). 

A suitable comparison sample for the SWELLS dataset 
is its parent sample, the SDSS-DR7 database itself We 
adopt the results of the recent structural analysis performed 
by L. Simard and collaborators (2011, in preparation). These 
authors flt two-dimensional models to the SDSS images, 
using the same parametrisation for bulge and disk as was 
adopted in this study. The range in redshift and stellar mass 
of their sample is very similar to ours, and the parent sam- 
ple is the same. To minimise differences with our sample, 
we construct the size mass relation for the SDSS sample 
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limited to h/a <0.6. Following iDutton et alj (|2010l ). we con- 
struct the size-mass relationship separately for the bulge, 
disk, and total size, adopting the semi-major ax i s size for 
all three components. In contrast to lDutton et alJ (|20ld ) we 
do not apply any selection based on stellar mass or spectral 
type. To minimise the impact of measurement biases we per- 
form the comparison by also using sizes and stellar masses 
estimated using SDSS images for our sample of lensefl The 
size-mass relation for SDSS and for our deflector galaxies is 
shown in Figure [6] The top panel shows the relation between 
the total half light radius and total stellar mass. The middle 
panel shows the relation between the half light radius of the 
disk and the total stellar mass. The bottom panel shows the 
relation between the half light radius of the bulge and the 
total stellar mass. As a population, the deflector galaxies are 
consistent with being drawn from the SDSS parent popula- 
tion. We quantify this statement by comparing the vertical 
offset of each data point with the mean and standard de- 
viation of the correlation, and computing a goodness of 
fit statisticlfl. We find ~ 26, 29, 9 respectively for total, 
disk, and bulge, with 27 degrees of freedom for total and 
bulge and 26 for disk (one galaxy was omitted because of 
an unphysically small size returned by the SDSS catalogue, 
< 1 kpc). This shows that our selection of targets at fixed 
stellar mass is unbiased with respect to the parent SDSS 
population once the b/a < 0.6 cut is taken into account. 
Clearly however, the distribution along the size-mass rela- 
tion is very different than that of the parent sample, since 
our target list was comprised almost exclusively of systems 
with mass above lO^"'^ M0 by design. 

It is also interesting to consider whether the distribu- 
tion of confirmed lenses is consistent with that of the parent 
population, or whether the confirmation rate appears to be 
a function of, e.g., stellar mass density. By comparing the 
solid (grade A lenses) and open points (grade B, C, X) in 
Figure [HI it appears that the solid points tend to lie pref- 
erentially somewhat below the open points. However, the 
difference is not significant. Even if we restrict to only the A 
grade lenses, they still appear to be consistent with the size 
mass relation of the SDSS parent population, at fixed stellar 
mass (X^ = 18, 19, 7 for 16 de grees of freedom, respectively 
for total, disk, and bulge component). 
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Figure 6. The size-mass relation for SWELLS deflector galaxies. 
The solid and dashed lines represent the mean, and 68% intrinsic 
scatter, of the size distribution of SDSS galaxies at fixed mass 
with b/a < 0.6. The points represent the SWELLS galaxies: solid 
squares are secure lenses (A), open stars are probable and possi- 
ble lenses (B and C), while open squares are not lenses (X). Red 
symbols identify pre-Cycle 16s targets, while black symbols iden- 
tify Cycle 16s targets. The top panel shows the total semi-major 
axis efi'ective radius as a function of total stellar mass. The mid- 
dle panel shows the semi-major axis efli'ectivc radius of the disk 
component as a function of total stellar mass. The bottom panel 
shows the semi-major axis efliective radius of the bulge compo- 
nent as a function of total stellar mass. We use the symbol R50 
to identify the semi-major axis effective radius as measured from 
SDSS images, to distinguish it clearly from the circularized half- 
light radius measured from HST images and given in Table [2] 



The SDSS sizes and masses are in remarkably good agreement 
to the HST determined ones, considering the difference in reso- 
lution. Nevertheless, the small differences are significant enough 
that they should be taken into account. Specifically, there is ex- 
cellent agreement between the SDSS and HST determined stellar 
mass values, with no significant offset (0. 04ifc0.03 dex), consistent 
with the findings of I Aueer etall | |2009|'| . In contrast, the HST- 
determined sizes are on average 0.10it0.04 (disk) and 0.17it0.06 
dex (bulge) smaller than those determined from SDSS images, 
with larger difl'erences for the smaller objects. These difl'erences 
are likely due to the higher spatial resolution of HST and/or a 
longer wavelength range used for our HST fits. 
* Errors on size can be neglected since they are much smaller 
than the intrinsic scatter of the correlation. 



6 FINDING SPIRAL LENS GALAXIES: 
LESSONS LEARNT AND FUTURE 
PROSPECTS 

Finding spiral lens galaxies is more difficult than finding 
early- type lens galaxies for two reasons: 1) spiral galaxies 
tend to be less massive and less dense than early-type gala:x- 
ies, and they therefore produce smaller image separations; 2) 
their complex surface brightness distribution and the pres- 
ence of dust makes it hard to identify multiple images if they 
are present. 

For these reasons a crucial element for the success of 
this search was the availability of multi-colour high resolu- 
tion imaging data. However, for any given candidate it is not 
clear a priori which band will be most useful. In general, the 
SLACS source s are by selection st ar-forming and therefore 
blue (see, e.g.. lNewton et ahlboill '). and stand out the most 
in the F450W filter. A good example of this phenomenon is 
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J1117+4704, where the quadruply imaged source near the 
centre gradually disappears moving to longer wavelengths. 
In some cases, however, the presence of dust makes the in- 
frared data crucial for studying the multiple images. An ex- 
ample of this phenomenon is J1703-I-2451, where the lensed 
source is barely visible in the optical images, while it is obvi- 
ous in the infrared one. Interestingly, two of the lenses were 
confirmed purely on the basis of K' AO imaging. It appears 
that multiband photometry is an essential ingredient for fu- 
ture studies. At optical wavelengths HST is the only option 
for the time being. With current technology AO is only fea- 
sible for ~l/3 of high latitude targets, and therefore a com- 
bination of AO (when possible) and HST (when AO is not 
possible) observations seems to be the most cost-effective 
strategy for the infrared. 

The other crucial ingredient is angular resolution. This 
is useful especially for the detection of compact sources or 
central counter-images. Limited resolution is one of the rea- 
sons why we do not confirm any lens with stellar Einstein ra- 
dius significantly smaller than one arcsecond. A good exam- 
ple of this phenomenon is J0841+3824, where the counter- 
image is clearly visible in /-band, but not at the lower resolu- 
tion i/-band, even though this may in part be due to colour 
variations. The WFPC2 images obtained in Cycle 16s signif- 
icantly under-sample the native HST PSF. Better-sampled 
images with ACS or WFC3-UVIS may result in an im- 
proved success rate. Higher resolution red/infrared images 
from Next Generation Advanced Optics or from the James 
Webb Space Telescope might combine the advantages of dust 
removal and angular resolution, even though colour contrast 
may still make blue/UV images competitive. Alternatively 
one may restrict the search to the systems with larger stellar 
Einstein radii to maximize success. However, high-resolution 
multiband images are needed to estimate this quantity ac- 
curately and therefore it does not appear a viable strategy 
for a large sample. Furthermore, one has to be mindful of 
the astrophysical implications of such a selection against ob- 
jects where dark matter contributes significantly to the cen- 
tral mass density or where the IMF normalization is higher 
than normally assumed. Not to mention that finding lenses 
with small Einstein Radii is astrophysically interesting in its 
own right, in order to probe the structure of the lower mass 
galaxies and those with smaller bulges. 

In terms of future prospects it seems that another way 
to improve the success rate would be to target emission lines 
in the sources at high resolution. This would be an efficient 
way to disentangle light from the source and defiector, and 
distinguish without any doubt multiple images from pecu- 
liar features in the surface brightness of the deflector. The 
rich fllter complement of WFC3 aboard HST makes it the 
most promising instrument for this purpose. Targeting in- 
frared lines might mitigate the effects of foreground extinc- 
tion, although with the present targets only Paschen a is 
accessible with decent strehl and the line is typically too 
faint for present day integral field spectrographs, even on 
large 8- 10m telescopes. A tunable filter or custom-made fil- 
ters will be necessary to pursue this goal. Alternatively, one 
could target lines in the radio mapping CO transitions or HI. 
ALMA and SKA will have the resolution and sensitivity to 
measure emission line maps of the source and lens, and will 
thus also offer the ability to measure kinematics of the lens 
and of the source. Molecular gas tends to trace the stars. 



so any source that is visibly lensed should also be lensed in 
CO. Neutral gas tends to be more extended (typically by a 
factor of 2 or more) than stars and molecular gas so this will 
increase the chances of strong lensing. 



7 SUMMARY 

The goal of the SWELLS (Sloan WFC Edge-on Late-type 
Lens Survey) is to study the distribution of luminous and 
dark matter in disk galaxies. With this goal in mind we 
have assembled an unprecedented sample of 27 candidate 
gravitational lens systems, where the defiector is a high in- 
clination disk galaxy. The sample is comprised of 8 systems 
previously discovered by the SLACS survey for which we 
have obtained follow-up multiband high resolution imaging, 
and 19 new targets. The new targets were selected from the 
SDSS database to have multiple redshifts in the SDSS fiber, 
as well as elongated morphology. We present multi-colour 
imaging of the sample, as well as structural parameters de- 
rived from this data for the whole sample. We then study the 
properties of the deflectors and compare them to a sample of 
non-lens galaxies selected from SDSS with the same criteria 
as our lens targets. Our main results can be summarised as 
follows: 

(i) Of the 19 new targets, 8 are classifled as A-grade (se- 
cure) lenses, 1 is classified as B-grade (probable), 6 are 
classified as C-grade (possible), and 4 are classified as X- 
grade (not lenses). Our global success rate in Cycle 16s was 
42 ± 13% for A-grade lenses. The total number of A-grade 
lenses in the SWELLS sample is 16. 

(ii) All the deflector galaxies have prominent disks ob- 
served at high inclination, by design. However, the sample 
spans a broad range of morphological and spectral types. It 
includes early-type spirals and lenticulars as well as later- 
type spirals, galaxies with strong Ha emission and galaxies 
with no detectable Ha emission. The stellar bulge to disk 
ratio of confirmed lenses range between 0.22 and 0.85, with 
an average of 0.55. 

(iii) All confirmed lenses have stellar masses above 10^°'^ 
M0 and stellar Einstein radii above 0'.'94. None of the low 
stellar mass density targets is confirmed to be a strong lens, 
consistent with the idea that stellar mass (as opposed to 
dark matter) dominates the mass density profile at the scales 
of a few kpc. 

(iv) The SWELLS targets and confirmed lenses follow the 
same size-stellar mass relation as a comparison sample of 
high-elongation {b/a < 0.6) non-lens galaxies selected from 
the SDSS survey. This implies that, at given stellar mass, 
SWELLS lenses are consistent with being representative of 
the overall population of high-elongation galaxies as identi- 
fied by SDSS. 

We conclude our study by summarising some of the 
lessons learnt during this search. Key elements for a suc- 
cessful disk lens search appear to be high resolution images 
spanning from the blue to the near infrared, to exploit colour 
contrast as well as low dust extinction. The main remaining 
limitation is our ability to disentangle light from the deflec- 
tor and source. Future studies should increase their success 
rate by targeting emission lines in the source at high resolu- 
tion, via narrow band imaging or integral field spectroscopy. 



20 Treu et al. 



ACKNOWLEDGMENTS 

We thank Matteo Barnabe for useful discussions and sug- 
gestions. TT acknowledges support from the NSF through 
CAREER award NSF-0642621, and from the Packard Foun- 
dation through a Packard Research Fellowship. A AD ac- 
knowledges financial support from a CITA National Fel- 
lowship, from the National Science Foundation Science and 
Technology Center CfAO, managed by UC Santa Cruz un- 
der cooperative agreement No. AST-9876783. AAD and 
DCK were partially supported by NSF grant AST 08-08133, 
and by HST grants AR-10664.01-A, HST AR-10965.02-A, 
and HST GO-11206.02-A. PJM was given support by the 
TABASCO and Kavli Foundations, and the Royal Society, 
in the form of research fellowships. LVEK acknowledges the 
support by an NWO-VIDI programme subsidy (programme 
number 639.042.505). This research is supported by NASA 
through Hubble Space Telescope programs CO- 10587, GO- 
11202, and CO-11978, and in part by the National Science 
Foundation under Grant No. PHY99-07949. and is based 
on observations made with the NASA/ESA Hubble Space 
Telescope and obtained at the Space Telescope Science In- 
stitute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 
5-26555, and at the W.M. Keck Observatory, which is op- 
erated as a scientific partnership among the California In- 
stitute of Technology, the University of California and the 
National Aeronautics and Space Administration. The Obser- 
vatory was made possible by the generous financial support 
of the W.M. Keck Foundation. The authors wish to recog- 
nize and acknowledge the very significant cultural role and 
reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We are most 
fortunate to have the opportunity to conduct observations 
from this mountain. Funding for the SDSS and SDSS-II was 
provided by the Alfred P. Sloan Foundation, the Participat- 
ing Institutions, the National Science Foundation, the U.S. 
Department of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS was managed by the Astrophysical 
Research Consortium for the Participating Institutions. The 
SDSS Web Site is |¥ttp:/7www.sdss.org/| 



APPENDIX A: SPECTROSCOPIC 
PRE-SELECTION OF SWELLS CANDIDATES 

The parent sample for the new SWELLS lens search consists 
of all SDSS DR7 spectra with either good or marginal plate 
quality that are classified by the SDSS spectroscopic reduc- 
tion pipeline as galaxies. Noise estimates for all spectra were 
rescaled on a plate-by-plate basis to be consistent with the 
statistics of sky-subtracted sky spectra on that plate. 

The spectroscopic selection algorithm of the original 
SLACS sample and our new SWELLS sample both rely 
upon subtracting a model for the continuum of the fore- 
ground galaxy and scanning the residual spectrum for higher 
redshift emission features. For SLACS, these models were 
taken directly from the SDSS redshift pipeline, and consist 
of the best-fit combination of 4 PCA-derived eigenspectra. 
For the absorption-dominated systems typical of SLACS, 



these models provide essentially noise-limited subtraction 
of the foreground galaxy continua. However, for emission- 
line foreground spectra, the quality of continuum subtrac- 
tion with these models is much worse. Hence, we derived 
new continuum models by (1) generating a new PCA basis 
with 7 eigenspectra, and (2) masking all common emission 
features in the foreground galaxy spectra before fitting each 
spectrum with this basis. 

The new models were subtracted from the target spec- 
tra, and the residual spectra were scanned automatically 
for either (1) multiple emission lines at a single background 
redshift, or (2) a single emission line consistent with a 
marginally resolved [OH] 3727 doublet. Potential [OH] iden- 
tifications were all checked for the presence of veto lines that 
would indicate Hq, [OIII]5007, or H/3 as a more likely iden- 
tification. The set of spectra identified by the automated 
scan was visually inspected by one of us (ASB) and pruned 
of obvious false-positive detections related to night-sky sub- 
traction residuals, broad emission-line wing residuals, and 
spectra of generally bad data quality. Neighboring fibers 
in the multifiber spectrograph image were also checked to 
make sure that candidate [Oil] detections were not due to 
crosstalk from bright emission lines in adjacent spectra, or to 
small-scale auroral emission in the night sky. The final list of 
potential target spectra was taken as the set which survived 
all visual inspections, and for which either (1) each of three 
distinct emission lines at the same background redshift was 
detected at 4-sigma or greater significance, or (2) a single 
candidate background [OH] emission feature was detected 
at 7.5-sigma or greater significance. Candidate [OH] systems 
were furthermore ranked as "A" or "B" based upon the ap- 
pearance of doublet line structure, which can be marginally 
resolved at SDSS resolution. 

In total, out of almost a million spectra in the SDSS 
archive, our search algorithm found more than 200 new high- 
probability lens candidates, in addition to those already tar- 
geted by previous SLACS programs (a similar number). The 
targets for the SWELLS program were selected from this 
new pool of candidates as described in the main text. 
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